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SUMMARY 


A  method  is  given  for  determining  the  horizontal  tail  loads  in 
maneuvering  flight.  The  method  is  based  upon  the  assignment  of  a  load- 
factor  variation  with  time  and  the  determination  of  a  minimum  time  to 
reach  peak  load  factor.  The  tail  load  is  separated  into  various  com¬ 
ponents.  Examination  of  these  components  indicated  that  one  of  the 
components  was  so  small  that  it  could  be  neglected  for  most  conventional 
.airplanes,  thereby  reducing  to  a  minimum  the  number  of  aerodynamic 
parameters  needed  in  this  computation  of  tail  loads. 

In  order  to  illustrate  the  method,  as  well  as  to  show  the  effect 
of  the  main  variables,  a  number  of  examples  are  given. 

Some  discussion  is  given  regarding  the  determination  of  maximum 
tail  loads,  maximum  pitching  accelerations,  and  maximum  pitching  veloc¬ 
ities  obtainable. 


INTRODUCTION 


The  subject  of  maneuvering  tail  loads  has  received  considerable 
attention  both  experimentally  and  theoretically.  Theoretically^  methods 
and  solutions  have  been  derived  for  determining  the  horizontal  tail  load 
following  either  a  prescribed  elevator  motion  (references  1  to  3)  or  an 
assigned  load— factor  variation  (reference  4). 

The  first  approach  has  been  adopted  into  some  of  the  load  require¬ 
ments  where  the  type  of  elevator  movement  specified  consists  of  linear 
segments  whose  magnitudes  and  rates  of  -  movement  are  governed  by  the 
assignment  of  a  maximum  initial  elevator  movement  consistent  with  the 
pilot* s  strength.  The  rates  of  movement  and  the  time  the  elevator  is 
held  before  reversing  are  so  adjusted  that  the  design  load  factor  will 
not  be  exceeded. 
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Tlie  results  of  reference  5  show,  as  Is  to  be  expected,  that  only 
when  the  aerodynamic  force  coefficients  are  accurately  known  from  wind- 
tunnel  tests  can  good  agreement  he  obtained  between  measured  and 
calculated  tail  loads.  At  the  design  stage,  however,  only  general 
aerodynamic  and  geometric  quantities  are  available  and  some  of  the  more 
important  stability  parameters  are  not  known  accurately.  Thus,  the  work 
involved  in  the  solution  for  the  tail  load  following  a  given  elevator 
motion  is  not  considered  to  be  in  keeping  with  the  accuracy  of  the  results 
obtained.  Consequently,  there  appears  to  be  a  need  for  an  abbreviated 
design  method  of  confuting  tail  loads  which,  although  incorporating 
approximations,  will  nevertheless  be  based  on  the  theoretical  considera¬ 
tions  of  the  problem. 

If  the  load-factor  variation  with  time  is  specified  and  the 
corresponding  tail  load,  elevator  angles,  and  load  distributions  are 
subsequently  determined,  a  simpler  and  equally  rational  approach  to  the 
tail— load  problem  can  be  made.  Although  this  approach  has  been  used  to 
a  limited  degree  (reference  4),  several  shortcomings  have  limited  its 
use. 


The  purpose  of  this  paper  is  to  develop  further  the  load— factor  or 
inverse  approach  and  to  present  a  method  of  computing  horizontal  tail 
loads  which  is  comprehensive  and  generally  simple.  To  this  end,  (1) 
the  shape  of  the  load— factor  curve  and  the  minimum  time  required  to 
reach  the  peak  load  factor  have  been  determined  from  an  analysis  of 
pull-up  maneuvers  that  were  available,  (2)  the  minimum  time  required  to 
reach  the  peak  load  factor  has  been  determined  from  a  theoretical  analysis 
which  is  supported  in  some  measure  by  statistical  data  obtained  from  a 
number  of  flight  tests  with  airplanes  of  widely  varying  sizes,  and  (3) 
the  equations  relating  the  various  quantities  are  presented. 


SYMBOLS 

b  wing  span,  feet;  also  shape  factor  in  equation  (13) 

bt  tail  span,  feet 

c  chord,  feet 

c  mean  aerodynamic  wing  chord,  feet 

Cl  lift  coefficient  (L/qS) 

Cm  pitching-moment  coefficient  of  airplane  without 

horizontal  tail  (Mb/qS^) 
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pitching-moment  coefficient  of  isolated  horizontal- 
tail  surface 

acceleration  of  gravity,  feet  per  second  per  second 

2 

pitching  moment  of  inertia,  slug-feet 

radius  of  gyration  about  pitching  axis,  feet 

empirical  constant  denoting  ratio  of  damping  moment 

of  complete  airplane  to  damping  moment  of  tail  -alone 

lift,  pounds 

local  lift  at  any  spanwise  station 
airplane  mass,  slugs  (W/g) 
moment,  foot-pounds 
airplane  load  factor  at  any  instant 
mn.-yi  rmim  increment  in  load  factor 
dynamic  pressure,  pounds  per  square  foot 


wing  area,  square  feet 
horizontal— tail  area,  square  feet 
time,  seconds 

time  to  reach  peak  of  elevator  deflection,  seconds 
airplane  true  velocity,  feet  per  second 
airplane  weight,  pounds 

length  from  center  of  gravity  of  airplane  to  aerodynamic 
center  of  tail  (positive  for  conventional  airplanes), 
feet 


nondimens ional  spanwise  dimension 
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constants  occurring  in  equations  (13)*  (23)*  (26),  and. 
(30) 

constants  occurring  in  "basic  differential  equation 
(see  equation  (3)  and  table  I) 

time  to  reach  peak  load  factor,  seconds 

mass  density  of  air,  slugs  per  cubic  foot 

tail  efficiency  factor  (q-fc/q)’ 

wing  angle  of  attack,  radians 

average  angle  of  attack  of  horizontal  stabilizer, 
radians 

tail  angle  of  attack,  radians 
angle  of  sideslip,  degrees 
flight-path  angle,  radians 
attitude  angle,  radians  (a  +  7) 
elevator  angle,  radians 

downwash  angle,  radians 


i^  tail  setting,  radians 

The  notations  co  and  0,  <x  and  0,  and  so  forth,  denote  single 
and  double  differentiations  with  respect  to  t. 

Subscripts; 

o  initial  or  selected  value 


t 

max 

*0 

geo 


tail 

maximum  value 
zero  lift 
geometric 


c 


camber 
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METHODS 

Method  of  Determining  the  Dynamic  Tail  Load 


Basic  equations  of  motion.—  The  single  differential  equations  for 
the  longitudinal  motion  of  an  airplane  for  any  elevator  deflection  (see 
method  given  in  reference  2)  may  "be  written  as 


rnyV 


Ax  qS  — 


i^qS^  A)  =  0 


.(1) 


dC_  o2  d-Cjn,  S+2 3  2  •• 

—  toq  -  tSL^  +  V  ^  if>  -  8-0 


(2) 


Equations  (l)  and  (2)  represent  summations  of  forces  perpendicular 
to  the  relative  wind  and  of  moments  about  the  center  of  gravity.  (See 
fig.  1  for  direction  of  positive  quantities.)  Implicit  in  these  equa¬ 
tions  are  the  following  assumptions: 

(1)  In  the  interval  between  the  start  of  the  maneuver  and  the 
attainment  of  maximum  loads,  the  flight— path  angle  does  not  change 
materially;  therefore,  the  change  in  load  factor  due  to  flight-path 
change  is  small. 

(2)  At  the  Mach  number  for  which  computations  are  made,  the  aero¬ 
dynamic  derivatives  are  linear  with  angle  of  attack  and  elevator  angle. 

(3)  The  variation  of  speed  during  the  maneuver  may  be  neglected. 

(k)  Unsteady  lift  effects  may  be  neglected. 

By  use  of  the  relations  9  =  y  +  <x,  0=7+  a,  and  6  =  7  +  a, 

equations  (1)  and  (2)  are  reducible  to  the  equivalent  second-order 
differential  equation 


ct  +  K-jO,  + 


K2  Ax,  =  K3  AS 


(3) 


where  Kq,  Eg*  and  K3  are  constants  for  a  given  set  of  conditions  (see 
table  I). 
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In  equations  (l)  and  (2)j  y9  Q9  AS>  and  AL^  will^  in  a 

given  maneuver 3  vary  with  time.  Using  the  relations  "between  y  y 

and  their  derivatives  permits  equation  (H)  to  he  rewritten  as  follows  to 
give  the  increment  In  tail  load: 


at,. 

t  da  hx^. 


o..  2“* 

mky  a.  inky  y 


_  -  b 

dS  btx^. 


In  a  -still  shorter  form,  equation  (4)  may  be  written  as 


AL+  =  AL.  +  AL,  +  AL.  .  +  AL . 
t  ta  ta  t7  tc 


Equations  (4)  and  (5)  show  that  the  tail-load  increment  (the  increment 

above  the  steady-flight  datum  value)  at  any  time  is  composed  of  four 

parts:  AL.  ,  associated  with  the  angle-of-attack  change;  AL,  ,  asso- 
ta  a 

ciated  with  angular  acceleration  about  the  flight  path;  AL.^,  associated 

with  angular  acceleration  of  the  flight  pa,th;  and  AL^.  ,  required  to 
~  c 

condensate  for  the  moment  introduced  by  change  in  camber  of  the  horizontal- 

tail  surface.  The  load  ALt  is  generally  small  but  in  some  extreme 

configurations  may  amount  to  10  percent  of  the  total  increment  and  thus 
for  the  present  it  is  retained  in  the  development. 

If  the  load-factor— increment  variation  with  time  An  is  known, 
then  by  the  usual  definition 


so  that 


Aa  = 


An  W/S 


h  W/S 

dOr 
- Ll  n 


n  W/S 
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The  following  relation  also  exists  "between  An  and  7: 


An  g  =  yY 


so  that 


7  = 


ng 


(7) 

(8) 


When  equations  (6)  to  (8)  are  substituted  into  equations  (4)  and  (5), 
the  four  tail— load  components  then  become 


ALt  =  Sr1  P-  An 

ta  dC^  fcxt 


W2ky2 

ALa.  =  -  - i—  n 

%  dC^ 


gSqxt 


(9a) 


(9h) 


da 


AL. 


n 


(9c) 


AL. 


d8 


nt*! 


A5 


(9d) 


Thus,  if  the  variation  of  the  load  factor  with  time  An  and  the 
geometric  and  aerodynamic  characteristics  of  the  airplane  were  known, 
the  first  three  components  of  the  tail  load  could  he  found  immediately. 

The  magnitude  of  the  fourth  component,  that  due  to  horizontal— tail  camber, 
would  follow  from  equation  (3)  in  which  the  elevator  angle  is  seen  to  he 


AS  = 


«  * 

a 


.  *2 
a  +  —  Ax 

K3 


(10) 


Substitution  into  equation  (10)  of  the  values  of  and  d 

from  equation  (6)  yields  the  value  of  the  elevator  angle  at  any  instant 
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A6  = 


w/s' 


E3  da 


(“ 


+  K±ii 


(id 


so  that,  finally,  the  fourth  eonponent  is  given  as 


AL 


4C“t 


dS 


\ 


w/s 


Vt  ^ 

da 


K- 


(n  +  Kxn  +  K2  Anj 


(12) 


The  procedure  outlined  shows  that  the  tail— load  magnitude  can  he 
determined  if  the  load-factor  variation  is  known. 

Types  of  load— factor  variation.—  The  relation  between  the  tail  load, 
•the  geometric  and  aerodynamic  characteristics,  and  the  load  factor  having 
been  established,  it  is  desirable  to  establish  a  load-factor  variation 
which  is  reasonable  as  well  as  critical  insofar  as  loads  are  concerned. 

The  maximum  value  of  load  factor  is  usually  specified;  however,  there 
are  many  possible  variations  for  the  shape.  Regardless  of  the  details 
of  shape,  the  load  factor  may  be  considered  to  rise  smoothly  and  con¬ 
tinuously  to  a  maximum,  the  rate  of  rise  depending  upon  several  variables. 
Beyond  the  maximum  value  of  the  load  factor  the  return  to  initial  condi¬ 
tions  can,  at  the  will  of  the  pilot,  be  either  gradual  or  rapid. 

Experiments  as  well  as  theoretical  studies  have  already  indicated 
that  the  maneuver  that  combines  maximum  angular  and  linear  accelerations 
causes  critical  loads  in  both  the  wing  and  tail.  One  such  maneuver  occurs 
when  the  maximum  load  factor  is  reached  as  rapidly  as  possible  by  using 
an  initial  elevator  movement  which  is  greater  than  that  required  to  reach 
a  given  steady-trim  value  of  the  load  factor.  This  initial  elevator  move¬ 
ment  is  followed  by  a  rapid  checking  of  the  maneuver  either  by  returning 
the  elevator  quickly  to  neutral  or  by  reversing  the  controls. 

The  shape  of  the  load-factor  curve  for  such  a  maneuver  may  be 
expressed  approximately  by  several  analytic  functions,  one  of  which  is 


An  =  at^e-0^ 


(13) 


By  way  of  illustration,  figure  2  shows  details  of  the  shape  of  the  load- 
factor  curve  obtained  with  the  use  of  equation  (13)  for  which  the  constants 
have  been  adjusted  so  that  an  8g  peak  is  reached  in  1  second.  By  further 
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adjustment  of  the  constants  the  load  factor  can,  within  certain  limits, 
"be  made  to  rise  to  any  specified  peak  and  to  diminish  in  any  prescribed 
manner. 

Because  the  positive  slopes  obtained  from  equation  (13)  are  always 
greater  than  the  negative  slopes,  the  positive  angular  accelerations  are 
greater  than  the  negative  ones.  In  general,  this  condition  is  true  for 
most  high  g  critical  maneuvers  performed  by  most  classes  of  airplanes, 
but  maneuvers  may  occasionally  be  performed  for  which  the  reverse  may  be 
true,  particularly  for  small  airplanes. 

Determination  of  constants.-  From  equations  (9),  (11),  and  (12) 
the  required  quantities  relating  to  load  factor  are  seen  to  be  An,  n, 
and  n.  Since  the  increment  An  is  to  be  given  by 


An  =  at^e~ 


then  at  maximum  load  factor 


"C) 

Thus  t  =  ^  at  maximum  load  factor.  Let  N  =  An  Then 


so  that 


(13) 

(1*0 

(15) 

(16) 

(17) 


Equation  (17)  is  in  nondimensional  form  where  X  is  the  time  to  reach 
the  peak  load  factor  and  b  is  a  constant. 
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When  equation  ( 17)  is  differentiated,  the  first  and  second  deriva- 
tives  'become 


and 


2£(i  _  1)  -  &  + 1 
t2\  hj  t 


iil2  _  An  ^2 
S  JS 


(18) 


(19) 


In  equations  (17)  to  (19)  the  quantities  N,  X,  and  b  are  now 
required  in  order  to  determine  the  variation  of  An,  n,  and  n.  The 
value  of  N  is  immediately  available  from  the  required  maneuver  load 
factor,  whereas  the  time  to  reach  the  peak  load  factor  X  can  he  ob¬ 
tained  from  examination  of  available  records  or  by  specification.  The 
constant  b,  as  may  be  seen  from  equation  (17),  can  best  be  described 
as  a  "shape"  factor  and  has  no  particular  physical  significance. 

The  values  of  X  and  b  should  be  associated  with  a  maneuver  which 
produces  maximum  tail  loads .  Therefore  the  time  X  to  reach  peak  load 
factor  should  be  the  minimum  possible  consistent  with  possible  pilot 
action  and  airplane  response.  The  shape  factor  b  should  also  be  con¬ 
sistent  with  both  of  these. 


In  connection  with  the  determination  of  the  minimum  time  to  reach 
peak  load  factor,  the  results  shown  in  figure  3  for  a  typical  airplane 
are  informative.  Figure  3(a)  shows  the  load-factor  variation  following 
several  abrupt  jump  elevator  movements.  The  load  factor  varies  with  the 
elevator  position,  but  the  time  to  reach  peak  load  factor  does  not. 
Figure  3(b)  shows  the  load— factor  variation  for  several  abrupt  hat- 
shape  elevator  impulses.  Again  the  load  factor  is  seen  to  vary  with 
the  amount  of  elevator  deflection  but  the  time  to  reach  the  peak  value 
remains  constant.  Although  the  time  to  reach  the  peak  load  factor 
shown  in  figure  3(b)  remains  constant,  it  is  seen  to  be  less  than  that 
shown  in  the  previous  case;  therefore,  an  inpulse  elevator  motion  pro¬ 
duces  a  smaller  value  of  X  than  the  jump  type. 

Because  of  inertia  and  elasticity  in  the  control  system,  the  pilot 
cannot  move  the  elevator  instantaneously  but  requires  some  finite  time 
tp  to  do  so.  A  possible  critical  type  of  elevator  impulse  thus  appears 

to  be  one  which  increases  linearly  to  maximum  and  decreases  at  the  same 
ra,te  to  zero.  In  order  to  determine  the  minimum  time  to  reach  peak  load 
factor  associated  with  such  a  variation,  the  equation  of  motion  (equa¬ 
tion  (3))  has  been  solved  for  the  triangular  elevator  impulse  for 
airplanes  of  various  static  stabilities  and  damping. 
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The  results  of  the  computations  are  given  In  figure  4  in  which  the 
minimum,  time  X  to  reach  peak  load  factor  is  plotted  against  the  time  t 
required  to  deflect  the  elevator. 

For  completeness  the  curves  of  figure  4  are  labeled  for  the  actual 
values  of  employed  in  the  computation  as  well  as  for  relative  values 

of  stability.  By  a  series  of  computations  the  damping  term  Kp  was 
found,  as  was  to  be  ejected,  to  have  only  a  secondary  effect  on  A,. 

The  curves  apply  to  an  average  value  of  the  damping  constant.  The  upper 
curve,  labeled  "low  stability,1'  should  be  associated  with  rearward 
center-of-gravity  positions  (that  is,  low  static  margin)  in  combination 
with  one  or  both  of  the  following:  low  dynamic  pressure  or  heavy  air¬ 
planes.  The  lower  curve,  labeled  "high  stability,"  would  be  associated 
with  forward  center-of-gravity  positions  in  combination  with  one  or  both 
of  the  following:  high  dynamic  pressure  or  light  airplanes.  It  is  seen 
that  X  increases  almost  linearly  with  tp  and  also  increases  when  the 
restoring  forces  are  reduced,  that  is,  when  the  stability  is  reduced. 

A  preliminary  value  of  the  shape  factor  b  (required  in  equa¬ 
tions  (17)  to  (19))  was  initially  determined  from  flight  records  of 
typical  impulse  maneuvers  by  fitting  curves  of  the  type  given  by  equa¬ 
tion  (13)  through  several  points  of  the  actual  time  histories  and 
determining  the  constants.  The  results  of  this  first  step  were  then 
modified  by  the  results  of  the  same  computations  which  had  been  made  to 
determine  X,  and  the  variation  of  b  with  tp  given  in  figure  5  was 
obtained.  Since  the  b  factor  is  not  found  to  be  critical,  an  average 
value  of  5.O  is  suggested,  although  as  a  refinement  the  values  from 
figure  5  may  be  used. 

The  question  of  the  value  of  tp  to  use  is  one  which  must  be 
solved  either  from  experience  or  from  a  knowledge  of  the  characteristics 
of  the  controls  and  the  control  system.  For  conventional  airplanes 
having  the  usual  amounts  of  boost  and  no  rate  restrictors,  the  following 
values  of  tp  are  suggested  as  representative: 


Fighters  or  small  civil  airplanes  with  weight  limit  from 

about  5OO  to  12,000  pounds,  seconds . .  0.20 

Two— engine  airplanes  with  weight  limit  from  25,000 

to  45,000  pounds,  seconds . .  0.25 

Four— engine  airplanes  with  weight  limit  from  50,000 

to  80,000  pounds,  seconds . . .  0.30 

Airplanes  with  weight  limit  above  100,000  pounds,  seconds  .  .  .'  0.40 
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The  minimum  time  X  given  in  figure  4  was  actually  established, 
separately  from  the  adopted,  load— factor  variation;  therefore,  in  applying 
the  inverse  method,  the  derived  elevator  impulse  would  not  be  expected 
to  agree  in  detail  with  the  "tent"  type  impulse  used  in  the  derivation. 

The  first  three  tail— load  components  can  now  be  computed  by  insert¬ 
ing  the  values  of  An,  h,  and  n  from  equations  (17)  to  (19)  into 
equation  (9)  and  using  appropriate  values  of  X  from  figure  4.  In ^ 
order  to  facilitate  this  computation,  curves  of  An/lf,  hx/N,  and  nXr/  N 
plotted  against  t/x  are  given  in  figure  6  for  the  suggested  value  of 
b  =  5.  Actually  to  apply  the  results  of  figure  6  it  is  convenient  to 
find  first  the  components  AL.  ,  and  so  forth,  in  terms  of  the  nondi— 

a 

mensional  time  t/x  and  then  to  convert  to  time  t  in  seconds.  In 
order  either  to  compute  the  fourth  component  or  to  obtain  the  elevator 
angles  for  use  in  chord  loading,  the  constants  Kp,  Kg,  and  K3  of 
equation  (3)  must  also  be  known. 

Thus,  in  terms  of  t/x  and  the  ordinates  of  figure  6,  the  various 
tail— load  components  are 


AL+  =  — ^  ^^.fordinate  of  fig.  6(a)) 
d(X  bx+\  > 


(20a) 


K^Ordinate  of  fig.  6(b))  . 

- - - -  +  Kg  (Ordinate  of  fig.  6(a)JJ(20d) 

The  constants  Kj_,  Kg,  and  K^  defined  in  table  I  are  the  same  as  those 
given  in  reference  2,  except  for  changed  signs  caused  by  specifying  xt 
as  positive. 
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The  conversion  to  time  t  is  made  by  multiplying  values  of  the 
base  scale  t/X,  by  \. 

Sample  calculations  for  incremental  tail  loads,—  The  results  of 
several  examples  are  given  to  illustrate  not  only  the  method  but  also 
the  effect  of  each  of  a  number  of  variables  on  the  incremental  tail  load 
of  a  typical  fighter  airplane,  the  geometric  and  aerodynamic  character¬ 
istics  of  which  are  given  in  the  following  tables.  In  order  to  illustrate 
the  effect  of  static  stability,  results  have  been  computed  for  three 
center-of-gravity  positions  with  the  assumption  that  an  8g  recovery  is 
made  at  19*100  feet  from  a  vertical  dive  at  an  equivalent  airspeed  of 
400  miles  per  hour.  In  order  to  illustrate  the  effect  of  the  time  of 
the  elevator  impulse  on  the  tail  load,  computations  were  carried  out  at 
one  of  the  center— of— gravity  positions  for  several  values  of  t^_.  The 

cases  considered  and  the  airplane  characteristics  follow: 


GEOMETRIC  CHARACTERISTICS 


Gross  wing  area,  S,  square  feet .  300 

Gross  horizontal— tail  area,  S-^,  square  feet . .  .  60 

Airplane  weight,  W,  pounds . . . 12,000 

Wing  span,  b,  feet .  4l 

Tail  span,  b^,  feet .  1 6 

Radius  of  gyration,  ky*  feet . ...... .  6.4 

Distance  from  aerodynamic  center  of  airplane  less  tail  to 
aerodynamic  center  of  tail,  x^_,  feet: 

Center  of  gravity,  29  percent  M.A.C .  20.0 

Center  of  gravity,  24  percent  M.A.C .  20.3 

Center  of  gravity  at  aerodynamic  center  ^  .  .  . .  21.0 


AERODYNAMIC  CHARACTERISTICS 


Slope  of  airplane  lift  curve,  dC-^/doc,  radians 
Slope  of  tail  lift  curve,  dC^/da^,  radians  .  . 

Downwash  factor,  de/da...  . 

Tail  efficiency  factor  ^yq^,  . 

Empirical  airplane  damping  factor,  K  . 

Elevator  effectiveness  factor,  dCj^/db,  radians 

Rate  of  change  of  tail  moment  with  camber  due  to 
elevator  angle,  dCm  /d6,  radian  •  ••••*. 


4.87 

3.15 

0.54 

1.00 

1.1 

1.89 

-0.57 


Ik 
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Rate  of  change  of  moment  coefficient  with  angle  of  attack 


for  airplane  less  tail,  dC^da,  radians: 

Center  of  gravity,  29  percent  M.A.C .  0.625 

Center  of  gravity,  2k  percent  M.A.C . .  0.1*03 

Center  of  gravity  at  aerodynamic  center . .  .  ...  0.000 


The  specified  conditions  for  the  sample  computations  are  given  in 
table  II.  The  computed  results  for  tail  components  are  given  in  fig¬ 
ures  7  and  8.  Figure  7  gives  results  for  varying  the  center  of  gravity 
and  figure  8  gives  similar  results  for  varying  t-]_.  The  tail— load 
components  are  computed  from  equations  (20)  and  the  derived  elevator 
angles  from  equation  (ll).  If  the  increment  in  tail  load  due  to  camber 
and  the  incremental  elevator  angle  are  not  required,  the  K  values  need 
not  be  computed  and  the  computations  are  considerably  shortened.  Fig¬ 
ures  7  and  8  show  that  a  maximum  error  of  only  about  4  percent  is 
introduced  by  this  omission. 


Method  of  Determining  the  Total  Tail  Load 


The  initial  or  steady— flight  tail  load  and  elevator  angles  to  which 
the  computed  incremental  values  are  to  be  added  must  also  be  determined. 
In  steady  flight,  the  horizontal  tail  furnishes  the  moment  required  to 
balance  the  moments  from  all  other  parts  of  the  airplane  so  that  the 
initial  load  may  be  written  as 


CmoqSc 

+ 


dCm  WS_ 
bxt 


cos  70 


(21) 


Thus  the  total  tail  load  at  any  time  in  a  maneuver  is  composed  of  the 
four  previously  mentioned  parts  plus  the  components  given  in  equa¬ 
tion  (21).  Only  the  first  term  of  equation  (21)  represents  a  new  type 
of  load  because  the  second  term  is  a  load  of  the  type  given  by  equa¬ 
tion  (9a)  or  equation  (20a)  and  its  effect  may  be  immediately  included 
in  the  computations  by  multiplying  the  ordinates  of  figure  6  by 
N  +  cos  7q  instead  of  by  N. 

The  initial  elevator  angle  required  to  balance  the  airplane  in 
steady  flight  varies  with  airplane  Cp  and  center— of— gravity  position 
so  that,  in  general,  60  must  be  obtained  from  wind-tunnel  data.  With¬ 
out  results  of  wind-tunnel  tests,  a  rough  rule  which  can  be  used  as  a 
guide  at  the  design  stage  in  determining  the  elevator  position  is  that 
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the  final  elevator  setting  will  he  so  adjusted  h y  repositioning  of  the 
stabilizer  setting  during  acceptance  tests  that  it  will  he  near  a  zero 
position  at  the  cruising  speed  and  at  the  most  prevalent  center— of— 
gravity  position. 


Method  of  Determining  Maximum  Yalues 

Maximum  tail  loads  and  angular  accelerations.—  The  method  outlined 
enables  a  point— by— point  evaluation  to  be  made  of  the  quantities  that 
determine  the  tail  load.  Such  detail  may  often  be  unnecessary  and  the 
procedure  may  be  shortened  by  evaluating  only  those  points  near  the  load 
peaks  or,  alternatively,  by  accepting  an  approximation  to  the  results. 

One  such  approximation  which  may  be  made  is  to  balance  the  airplane  at 
the  combinations  of  load  factor  and  angular  acceleration  which  would 
result  in  maximum  up  and  down  tail  loads. 

Figure  7  shows  that  the  maximum  down  tail  load  in  a  pull-up  occurs 
near  the  start  of  the  maneuver  and  before  appreciable  load  factor  is 
reached.  This  maximum  load  is  practically  coincident  with  the  negative 

maximum  in  the  L+.  tail— load  .component.  Since,  for  a  given  configu— 

a 

ration,  this  component  increases  as  the  center  of  gravity  is  moved  forward 
and  since  the  steady— flight  down  load  increases  with  speed,  the  maximum 
down  tail  load  in  a  pull-up  occurs  at  the  highest  design  speed  in  combin¬ 
ation  with  the  most  forward  center— of-gravity  position. 

Figures  7  and  8  show  that  at  the  time  of  the  maximum  down-tail— load 
increment  the  elevator  is  near  but  has  not  quite  reached  its  peak  position. 
Also  at  the  time  of  maximum  up— tail— load  increment  the  elevator  is  near 
its  zero  position,  although  it  may  be  on  either  side  of  this  position 
depending  upon  the  stability  and  the  time  tj_.  These  results  suggest 
that  the  maximum  down  load  for  the  elevator  and  the  hinge  brackets  would 
occur  with  the  airplane  center  of  gravity  well  forward  and  at  the  start 
of  the  maneuver .  The  maximum  load  for  the  stabilizer  is  likely  to  occur 
at  the  peak  load  factor. 

Figure  7  also  shows  that  the  up  tail  load  occurs  near  the  peak  of 
the  component  as  well  as  near  the  positive  maximum  peak  in  the 

component.  Since  the  •  component  increases  as  the  center  of 

gravity  is  moved  rearward  and  since  a  decrease  in  speed  generally  reduces 
the  initial  down  load,  the  maximum  up  tail  load  occurs  at  the  upper  left- 
hand  corner  of  the  Y— n  diagram  for  the  most  rearward  center-of-gravity 
position. 
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The  maximum,  tail  load  in  a  pull-up  maneuver  may  be  written  as 


dO. 


m  WS 


max 


(n  +  cos 


7o)~if 


(22) 


where  the  sum  of  the  second  and  third,  terms  is  to  "be  a  maximum  in  the 
maneuver.  From  the  previous  discussion  the  load— factor  increment  at 
maximum  down  load  is  nearly  zero  and  at  maximum  up  load  it  is^nearly 
equal  to  N  so  that  if  the  positive  and  negative  values  of  can 

he  determined,  a  relatively  simple  method  for  determining  maximum  loads 
is  available. 

Since  by  definition  0  =  a  +  7,  an  expression  for  angular  acceler¬ 
ation  can  be  derived  from  equations  (6)  and  (7)  and  written  in  the  form 


•-  w/s  .. 

9  =  tt —  n 
°L  4 


The  maximum  angular  acceleration  can  be  approximated  by 


9  .  B  -S£-  1- 

max  *  ^2 

da  ^ 


+  C 


vx 


(23) 


For  the  maximum  positive  pitching  acceleration,  B  is  the  maximum 
positive  ordinate  in  figure  6(c)  and  C  is  the  ordinate  of  figure  6(b) 
at  a  value  of  t/x  for  which  B  was  determined.  Thus,  B  is  6.5  and 

C  is  0*95  for  this  example. 

For  the  maximum  negative  pitching  acceleration,  B  is  the  maximum 
negative  ordinate  in  figure  6(c)  and  C  is  the  ordinate  of  figure  6(b) 
at  a  value  of  t/x  for  which  B  was  determined.  Thus,  B  is  —5.8  and 

C  is  0.80.  For  use  in  equation  (23)  the  values  of  X  for  the  maneuver 

are  available  from  figure  4  and  the  other  quantities  are  available  from 
the  conditions  of  the  problem.  The  maximum  loads  can  be  given  by  the 
following  equations : 
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For  maximum  up  tail  load  in  the  pull—' up: 


Jmax+  xt 


%  +  5s®.  (N  +  i.o)  +  hx(5JL2l&-  0i8a 

IOl  *>xt  ^  +  J-UJ  xt  MdCL  ,  V 

da  qA" 


(24a) 


For  maximum  down  tail  load  in  the  pull-up: 


■’max— 


S.^0  ,  jgm  WS 
Xj.  dCb  bx+ 


da  ^ 


(24h) 


M 


For  push-downs  to  limit  load  factor,  equations  (24a)  and  (24b)  still 
apply  with  changed  signs  for  W  and  changed  directions  for  and 

max+ 

L-l.  .  A  question  arises  as  to  whether  the  maximum  down  tail  load  at 
max— 

the  start  of  a  pull-up  with  forward  center— of-gravity  position  is  greater 
than  that  which  would  occur  when  pulling  up  from  a  negative  load— factor 
condition,  with  the  center  of  gravity  in  the  most  rearward  position. 

This  can  be  determined  only  by  computing  both  cases  and  seeing  which  is 
the  larger . 


Maximum  value  of  angular  velocity.—  The  maximum  value  of  the  pitching 
angular  velocity  in  the  pull-up  may  also  be  found  in  a  manner  similar  to 
that  used  to  obtain  the  maximum  angular  acceleration.  Since  9  =  a  +  y 
and  the  relations  involving  these  quantities  in  terms  of  load  factor  are 
given  by  equations  (6)  and  (7),  the  following  equation  may  be  written: 


•  n  W/S  An  g 

S  -  -J-  +  — 


(25) 


The  maximum  angular  velocity  may  be  approximated  by 


6  = 


X  dCL 


+  E 


Kg 

Y 


(26) 
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vhere  D  is  the  maximum  positive  ordinate  in  figure  6(h)  and  E  is  the 
ordinate  of  figure  6(a)  at  a  value  of  t/x,  for  which  D  was  determined. 
Thus  D,  for  this  example,  is  1-95  and-  E  is  0-4-8* 

Tn  the  steady  turn  or  pull-up  at  constant  g,  the  angular  velocity 
is  usually  given  hy  the  expression  9  =  1.0  Sp  The  difference  between 

the  factor  1.0  of  this  expression  and  the  factor  0.48  of  equation  (26) 
is  more  than  made  up  by  the  angle-of-attack  component  of  the  angular 
velocity. 


Approximate  Method  of  Determining  Load  Distribution 

g ymrnA-h-ri gal  i  oad  1  ng .  -  The  spanwise  distribution  of  the  total  load 
can  be  formulated  with  various  degrees  of  exactness.  If  information 
regarding  details  of  the  angle— of— attack  distribution  across  the  span 
were  available,  then  an  exact  solution  could  be  obtained  for  the  loading 
with  the  use  of  existing  lifting-surface  methods.  The  following  method 
may  be  used  as  a  first  approximation  to  the  solution. 

From  the  total  tail  load,  the  total  tail  lift  coefficient  can 

readily  be  found.  The  average  effective  angle  of  attack  a  of  the 
stabilizer  portion  is  given  in  the  definition 

-  J f1  tdJ*  *  L  °!5(S°  +  ^  I  dy*  (£7) 


where  only  a  Is  assumed  as  unknown  and  c7  and  c7  may  he  taken 

.  ■ a  o 

as  the  rates  of  change  of  section  lift  coefficient  with  <x  and 
respectively. 

Thus.,  for  constant  elevator  angle  across  the  span., 
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In  a  practical  case  both,  integrals  in  equation  (28)  need  he  evaluated 
only  once  for  a  given  configuration  and  Mach  number.  A  plot  of  & 
against  with  6  as  a  parameter  would  be  useful  in  further  compu- 

"fc 

tations.  With  a  known  as  a  function  of  and  5,  the  local  lift 

t 

at  any  spanwise  station  is  then  obtained  from  the  expression 


Z  =  c^qc  =  a  +  Cj6(80  +  AS)]  qc 


(29) 


IJhjgymmetrlcal  loading.—  Up  to  this  point  the  total  loads  have  been 
assumed  to  be  symmetrical  about  the  airplane  center  line,  whereas,  in 
reality,  the  load  may  have  an  asymmetrical  part.  The  sources  of  this 
dissymmetry  may  be  due  to  uneven  rigging,  differences  in  elasticity 
between  the  two  sides,  or  to  effects  of  slipstream,  rolling,  and  sideslip 
The  first  two  sources  are  usually  inadvertent  ones  while  the  last  two  are 
difficult  to  determine  without  either  wind-tunnel  tests  or  a  knowledge  of 
how  the  airplane  will  be  operated.  Present  design  rules  regarding  dis¬ 
symmetry  of  tail  load  are  concerned  more  with  providing  adequate  design 
conditions  for  the  rear  of  the  fuselage  than  with  recognizing  that  at 
the  maximum  critical  tail  load  some  dissymmetry  may  exist. 

Tests  in  the  Langley  full-scale  tunnel  (reference  6)  and  flight 
tests  (reference  7)  of  a  fighter-type  airplane,  as  well  as  unpublished 
flight  tests  of  another  fighter^-type  airplane,  indicate  that  the  tail- 
load  dissymmetry  varies  linearly  with  angle  of  sideslip  so  that  the 
difference  in  lift  coefficient  between  the  two  sides  of  the  tail  can  be 
given  as 


CL+  ~CL+  =  W  (30) 

Tttght  TLeft 

The  average  values  of  A  per  degree  found  for  the  two  fighter— type 
airplanes  are  approximately  0.01.  No  similar  values  are  available  for 
larger  airplanes  nor  for  tail  surfaces  having  appreciable  dihedral. 

In  maneuvers  of  the  type  considered  herein  it  is  doubtful  that 
angles  of  sideslip  larger  than  3°  would  be  developed  at  the  time  the 
maximum  tail  load  is  reached.  If  the  value  of  the  sideslip  angle  at 
the  time  of  marl  mum  tail  load  can  be  established,  equations  (27)  to  (29) 
are  easily  modified  to  include  this  effect,  provided  the  approximate 
value  of  A  is  known. 
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Chordwise  loading*  —  The  chordwise  distribution  can  "be  determined 
for  any  one  spanwise  station  in  either  of  two  ways.  One  way  for  design 
work  is  outlined  in  reference  8.  A  knowledge  of  the  airfoil  section 
and  the  quantities  contained  in  equation  (29)  suffices  for  this  deter¬ 
mination. 

If  pressure— distribution  data  are  available  for  a  similar  section 
with  flaps *  an  alternate  way  would  he  to  distribute  the  load  chordwise 
according  to  the  two-dimensional  pressure  diagrams  with  the  use  of  the 
computed  values  of  section  lift  coefficient  and  elevator  angle. 


DISCUSSION 


The  method  presented  is  another  approach  to  the  determination  of 
tail  loads.  From  the  results  given  in  figures  7  an<l  8/  it  can  be  seen 
that  the  camber  component  L^.  is  so  small  that  for  all  practical  cases 

it  may  be  omitted  with  considerable  simplification  In  the  computation  of 
tail  loads.  This  omission  reduces  to  a  m.lnlimim  the  number  of  aerodynamic 
parameters  needed  to  compute  the  tail  loads. 

It  is  possible *  in  the  application  of  the  present  method  with  the 
use  of  the  suggested  values  of  tj_*  that  the  derived  elevator  angles 

may  not  be  within  the  pilot* s  capabilities.  Since  it  must  be  assumed 
that  all  airplanes  *  to  be  satisfactory*  should  have  sufficient  control 
to  reach  their  design  load  boundaries*  such  an  occurrence  requires  only 
that  the  time  to  reach  elevator  peak  deflection  tq  be  increased  so  as 

to  reduce  the  elevator  angle.  The  results  of  figure  8*  in  which  the 
time  t]_  is  varied*  furnish  a  useful  guide  for  determining  the  in¬ 
crease  tq  that  might  be  required. 

If  sufficient  information  is  available*  it  is  recommended  that 
existing  lifting-surface  methods  be  used  in  determining  the  spanwise 
distribution  of  the  total  load;  however*  if  information  of  the  angle- 
of— attack  distribution  across  the  span  is  not  known*  the  method  presented 
may  be  used  as  a  first  approximation. 

Along  some  of  the  boundaries  of  the  V-n  diagram,  tail  buffeting  may 
occur.  Measurements  show  that  buffeting  usually  occurs  along  the  line 
of  maximum  lift  coefficient  and  again  along  a  high-speed  buffet  line 
which  is  associated  with  a  compress ibility  or  force  break  on  some  major 
part  of  the  airplane.  All  airplanes  are  subject  to  buffeting  at  the 
design  conditions  associated  with  the  lefi>-hand  comer  of  the  V-n 
diagram. 
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Only  high-speed  and/or  high-altitude  airplanes  are  capable  of  reaching 
the  other  boundary.  Measurements  show  that  the  oscillatory  buffeting 
loads  may  be  so  high  that  the  designer  should  at  least  be  cognizant  of 
them  at  the  design  stage. 

The  maximum  angular  acceleration  varies  inversely  "with  airspeed 
and  directly  with  the  load  factor*  with  the  contribution  due  to  acceler¬ 
ation  in  angle  of  attack  likely  to  be  more  important  than  the  angular 
acceleration  of  the  flight  path.  A  somewhat  similar  variation  is  indi¬ 
cated  for  the  maximum  angular  velocity  (equation  (26))  where  it  is  seen 
by  direct  substitution  that  the  part  due  to  angle  of  attack  is  likely 
to  be  larger  than  the  part  due  to  the  angular  velocity  of  the  flight 
path. 


CONCLUDING-  REMARKS 


A  simple  method  has  been  presented  for  determining  the  horizontal 
tail  loads  in  maneuvering  flight  with  the  use  of  a  prescribed  incremental 
load— factor  variation. 

The  incremental  tail  load  was  separated  into  four  components  repre¬ 
senting  a*  a*  y *  and  c.  The  camber  component  L-j.  is  so  small  that 

c 

for  most  conventional  airplanes  it  may  be  neglected*  thereby  reducing  to 
a  minimum  the  number  of  aerodynamic  parameters  needed  in  this  computation 
of  tail  loads. 

An  approximate  method  is  presented  for  predicting  maximum  angular 
accelerations  and  maximum  angular  velocities. 

The  method  indicates  that  maximum  tail  loads  in  a  pull— up  occur  at 
forward  centen-of-gravity  positions  and  early  in  the  maneuver.  The 
maximum  down  tail  loads  In  a  pull-up  occur  at  the  highest  design  speed 
in  combination  with  the  most  forward  centen-of— gravity  position.  The 
maximum  up  tail  load  occurs  at  the  upper  left-hand  comer  of  the  V— n 
diagram  for  the  most  rearward  center-of— gravity  positions. 
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TABLE  II 

SPECIFIED  COUDITIONS  OF  SAMPLE  PROBLEM 


Increment  in  load  factor  . . . . '  •  •  0 

Altitude,  feet . .  19,100 

Air  density,  slug  per  cubic  foot  .  . . . O.OOI3O6 


Case 

Gg 

(percent  M.A.  C. ) 

*i 

Ki 

*2 

e3 

X 

(fig.  b) 

1 

a  •  c  • 

0.2 

^.93 

30.4 

-33.4 

0.45 

2 

24 

.2 

4.72 

16.2 

-32.2 

•  5° 

3 

29 

.2 

4.61 

8.45 

-31.7 

.56 

24 

.4 

4.72 

16.2 

-32.2 

•77 

5 

24 

.6 

4.72 

16.2 

-32.2 

1.02 
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Figure  2.-  Variation  of  load-factor  increment.  An 
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Figure  5*-  Variation  of  shape  factor  b  with  t 
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Figure  6.-  Variation  of  incremental-load- factor  curves  with  time  ratio 
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Figure  7*-  Effect  of  center-of -gravity  position  on  incremental-tail-load 

components. 
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Figure  8.-  Effect  of 


T/me  r'Sec 

on  incremental -tail -load  components.  Center 
of  gravity,  0.24  c. 
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